Highest weight state description of the isotropic spin-1 chain 
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We introduce an overcomplete highest weight state basis as a calculational tool for the description 
of the isotropic spin-1 chain with bilinear exchange coupling Ji and biquadratic coupling J2. The 
ground state can be expressed exactly at the three special points in the phase diagram where the 
Hamiltonian corresponds to a sum of nearest neighbor total spin projection operators (Ji = > J2, 
Ji = — J2 < 0, and Ji = — J2/3 > 0). In particular, at the phase transition point Ji = — J2 < it 
is possible to exactly compute the ground states, excited states, expectation values, and correlation 
functions by using the new total spin basis. 



I. INTRODUCTION 

There has been a large interest in the isotropic 
one-dimensional spin-1 chain ever since Haldanc's 
prediction? 1 ^ that the excitation spectrum in integer spin 
Heisenberg chains should show a gap in strong contrast 
to the model with half integer spins. The general SU (2) 
invariant spin-1 chain model with nearest neighbor cou- 
pling is given by 

N 




H = ^ ( JiSi • S i+ i — J 2 (Sj • Sj + i) 2 ) 

i=l 
JV 

= J^ (coses, -S i+ i -sin9(Si ■ S m ) 2 ) , (1) 



where S,; are the spin-1 operators at the site i in 
a one-dimensional periodic lattice with N sites. Ex- 
act analytical solutions at special points were ob- 
tained by Affleck, Kennedy, Lieb and Tasaki (AKLT)^ 
Sutherland^ Kliimper^iii^ Barber— which supported 
Haldane's hypothesis and established an interesting 
phase-diagra m 10 ! 11 as shown in Fig. [TJ Experimen- 
tal results on quasi-one-dimensional spin-1 compounds 
such as "NENP" , CaNiCl 3 or AgVP 2 S 6 also confirmed 
the gap and the existence of effective spin-i spins near 
boundaries. 12 ' 13 The spin-1 chain was also one of the driv- 
ing forces in the development of the density matrix renor- 
malization group (DMRG) algorithm which in turn pro- 
vided very accurate estimates of the excitation spectrum 
and the correlation lengths at the Heisenberg point i 14 ' 15 
By changing the ratio of the Heisenberg coupling J\ 
and biquadratic exchange term J2, the system can be 
tuned through at least three antiferromagnetic regions 
and one ferromagnetic phase as shown in Fig. [1] The an- 
tiferromagnetic phase consists of three regions, Trimer, 
Haldane, and Dimcr, of which the last two are gapped. 
Because in most substances the biquadratic exchange 
term is much smaller compared to the bilinear term, 
the experimental realization in regions with dominant bi- 
quadratic exchange term was not possible for a long time. 
First experimental success was achieved with LiVGe2 02 
which appears to be well described by a large positive 
value of J2-— 

The similarity of the phase diagrams of the spin-1 chain 



FIG. 1: lhs; Phase diagram of the spin-1 chain: A: Pure 
Heisenberg chain; B: AKLT point [33; C: Sutherland 
model 51; D: Phase transition 19]; E: Phase transition 
mUL.? 1 : Exactly solvable 0011; G: SU(2) 2 integrable 
point 22 .23.24]. rhs: Phase diagram of the spin-^ chain: 
a: Heisenberg chain [25[; b: Majumdar-Ghosh point [13]; c: 
Critical point [26l ]: d: Phase transition; g: Two independent 
chains, St/(2)i x SU(2)i integrable; e: Phase transition [2?J. 
The points B, D, and E are "projection points" treated in 
this paper. 



compared to the spin- ^ chain with next nearest neighbor 
coupling is striking. In particular, the spin-^ chain also 
shows three antiferromagnetic regions, two of which are 
believed to be gapped and one ferromagnetic phase as 
shown in Fig. [T] It is known that the AKLT point in 
the spin-1 chain is in the same phase as the Majumdar- 
Ghosh pointi in the spin-^ chain, i.e. the two points B 
and b in Fig. Q] can be connected in a more general pa- 
rameter space £L We also see that the two gapped phases 
in both the spin-1 and the spin- ^ chain are separated by 
an integrable point with SU(2) 2 symmetry (points G, g). 
Moreover, in both cases the phase transition to the fer- 
romagnetic behavior occurs at points where the Hamil- 
tonian can be written in terms of total spin projection 
operators of neighboring spins (points D, d, E, e). There- 
fore, the main difference between the two phase diagrams 
is that in the spin-1 chain the isotropic point A happens 
to be in a gapped phase, while the isotropic point a falls 
in a gapless region in the spin-i chain. 

In this paper, we develop a highest weight state basis 
which can be used to calculate a number of exact prop- 
erties, especially at the "projection points" B, D and E, 
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where the Hamiltonian can be written as a sum over sin- 
glet, triplet, or quintet projection operators of two neigh- 
boring spins, respectively. In section HTl we introduce the 
general total spin basis, which is in spirit similar to the 
valence bond basis for spin-i systems where the total 
spin of pairs of spins in the system is specified. In section 
IIII1 we demonstrate how to calculate the ferromagnetic 
excitations with total spin s = N — 1 in the new basis 
as a simple illustration. At the singlet projection point 
D, the ferromagnetic excitations become dispersionless 
and it is possible to construct s = N — 2 excitations in 
agreement with earlier results^ as shown in section IIVI 
In section fVl we consider the triplet projection point E 
'At J i = — J 2 < (O = j) in a system with an even 
number of sites N. At this phase transition the antifcr- 
romagnctic s = and ferromagnetic s = N ground states 
are degenerate, but also ground states with any even to- 
tal spin s exist, leading to a large degeneracy. We are 
able to calculate the correlation functions exactly which 
remarkably do not decay along the chain, even for the an- 
tiferromagnetic ground state. Using the highest weight 
state basis, excited states can also be constructed exactly 
at this point. For completeness we also show how to ex- 
press the AKLT state at the quintet projection point B 
(G = — arctanl/3) in the new basis in section IVTl We 
conclude in section IVTll The appendix explains how op- 
erators are applied and scalar products are computed in 
the highest weight state basis. 

II. THE TOTAL SPIN BASIS FOR SPIN-1 
SYSTEMS 

In this section a basis set of highest weight states with 
given total spin s is introduced. Similar to the valence 
bond basis for spin-i systems, it is possible to specify 
pairs of spins (so-called valence bonds) which have defi- 
nite total spin. In spin-1 systems, pairs of spins can now 
be quintets (s = 2), triplets (s = 1) or singlets (s = 0). 
The most important difference to an ordinary S^-basis is 
that in this construction " bonds" are specified instead of 
local quantum numbers. 

We use a notation similar to the conventional 
one for the spin-i chain by Majumdar and Ghoshii 
{«i,«2}m=o = |s = 0,m = 0), [ii,i 2 ]m = \s = l,m) and 
(*i;*2)m = \s = 2, m), where s is the total spin of the 
pair of spins at sites i±, and m the S z eigenvalue. An 
important simplification for SU (2) invariant systems is 
that only highest weight states (hws) with m = s have 
to be considered as representatives of degenerate multi- 
plets, since it is possible to construct an entire degener- 
ate multiplet by applying the total spin lowering operator 



S~ = , which commutes with H. 

Accordingly, it is useful to define "bonds" as pairs of 
two spins, which are in a highest weight state of total 
spin 2, 1, or 

(ii, 12)2 = (ii,«2) = |++) 
[ixMi = [»i,*2] = (|+0>-|0+))/>/2 
{H,i 2 } = {i lji2 } = (|+_) + |_+)_|00))/V3.(2) 

Singlet and quintet bonds have even parity under ex- 
change of indices, while the triplet bonds are odd. Using 
this notation it is now possible to construct a highest 
weight state of the total system by specifying all bonds 
as follows 

Q t s 

\ips{h,-,iN)) = (h,h) ■■■(;■) [;■] {•,-}-{ijv-i,*Jv} , 

(3) 

assuming an even number of lattice sites TV. The number 
of quintet-, triplet- and singlet-bonds, Q, T and S defines 
the total spin s = 2Q + T. The wave-function is a func- 
tion of the pairwise different ordered indices i\, 1%, «at, 
which arc used to specify the bonds. The order of the 
quintet indices does not matter since all spins in quintet- 
bonds are in the state |+), so Eq. © can be simpli- 
fied by only specifying which spins belong to the quintets 
(ix, «2q)- Also the order within each singlet and triplet 
bracket is irrelevant up to a possible minus sign. To sim- 
plify the calculations it is sometimes useful to specify the 
bonds graphically as follows 



{hih}=/^\ [ii,i2]=r (ii,i 2 ) = •• , (4) 

where no bond connections are assigned to quintets. 

Clearly, states of the form (J3j) are indeed highest weight 
states of the system, since 

Stot IVO = m\ip a ) 

s? ot ivo = 8(s + 1) m , (5) 

where s = m = 2Q + T. Here, the first line is obvious 
from Eq. ([2]) and the second line follows since the total 
spin s has to be at least m, but can be at most the sum 
of the individual bond spins s = 2Q + T. 

States of the form ([3]) can now be used as basis states in 
order to express any hws of the system. For this purpose 
it is sufficient to only consider basis states with no triplets 
T = (for s = 2Q even) or one triplet pair T = 1 (for 
s = 2Q + 1 odd), since states with T > 1 can be expressed 
by the following linear combination 
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[l,2][3,4]tx (2,3){1,4} +(1,4){2,3}-(1,3){2,4}-(2,4){1,3} 

.r\. - . 



(6) 

In the hws basis defined in this way, it is possible to determine the action of any SU (2) invariant operator involving 
two spins, i.e. linear combinations of S; • S; + i and (Si • Sj+i) . For example for states with two singlet bonds we find 
the following relations 



or equivalently 



Si-S i+1 {j,i}{i + l,k} = {j,i + l}{i,k}- {j,k}{i,i + 1} 
(S t -S l+1 ) 2 {j,i}{i+l,k} = {j,i}{i + l,k} + {j,k}{i,i+l} , 



(7) 



J ii- 



1 k 



j % i 



1 k 



( Si .s i+l) 2 r\ r\ = r\ r\ 



j i i + 1 k 



J i 1 



1 k 




(8) 



Since the terms in the Hamiltonian (TT]) operate only on 
two spins at a time, it is always sufficient to consider 
clusters of two bonds irrespective of the length of the 
chain. A complete list how Si • Sj+i and (Si • Si + i) 2 
operate on the states can be found in appendix [X] 

From those relations (jXJ it is clear that the represen- 
tation of linear combinations in the hws basis in Eq. ([3]) 
is closed under the operation of any local SU(2) oper- 
ator, including permutations. Therefore, starting with 
any hws all other states in the corresponding sector of 
the Hilbcrt space can be represented as a linear combi- 
nation in the hws basis (J3j) . However, the linear combi- 
nations are not necessarily unique since the hws basis is 
overcomplete and not orthogonal. In particular, for the 
s = sector the hws basis is linearly independent up to 
N = 6 spin, while for N — 8 the basis becomes over- 
complete. This means that no relation between singlet 
states with 3 singlet bonds exists, while there are 14 (rel- 
atively complicated) relations involving 4 singlet bonds, 
which will not be discussed here. In general, there are 
iV !/(2( iV - s )/ 2 (i^)!s!) basis states of the form © for a 
given even spin s = 2Q. For comparison, if only the 
quantum number in the z-direction S z = in the ordi- 
nary S z basis is specified, the number of basis states is 

given by l+EiV+=i (jv- 2 iv^)!iv + !JV + ! which is a lar S er ba " 
sis than the s = hws basis up to N = 18. The situation 
for N = 4 spins is discussed in detail in appendix [B] as 
an example. 

It is also possible to determine scalar products in the 
hws basis using a straight-forward algorithm as presented 
in appendix [Cl In particular, the scalar product be- 
tween two total spin s = states can be obtained by 
the minimum number 7 of index exchanges Vi k j ; needed 
to transform the indices i\, ijy into an equivalent bond 



configuration corresponding to indices i[, ■■■ 1 i / N 

(ip s =o{h, ■■■,iNMs=o(i'i, -,i' N )) = ( 9 ) 

7 

(lpa=o[h,-,iN)\ V...V 1^8=0 («i,-,*jv)) = ^7 • 



Interestingly, this means that any basis state in the s = 
sector is never orthogonal to any other s = basis state. 



III. FERROMAGNETIC EXCITATIONS 

In order to demonstrate how to apply the new basis, we 
consider as an example excitations on the ferromagnetic 
state \F) = (1,2, ...,7V) which is always an eigenstate 
with energy E Q = N(Ji - J 2 ) = J/V(cos0 — sin 9). In 
order to construct a spin- wave excitation \k) with total 
spin s = N — 1 we can write 

\k) = e^[h,i2](i3,...,iN) • (10) 

Here, we require k ^ since the k = "excitation" of the 
ordinary spin-wave construction actually corresponds to 
a s = N multiplet state. We will now show explicitly that 
the state in Eq. PU|) is an eigenstate of the Hamiltonian 
(fT|) for any J\ and J 2 . In order to apply the operators 
Si • Sj+i and (Si ■ Si + i) 2 , we write states involving triplet 
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bonds with the sites i or i + 1 separately 

= e lk %i + \]{...) + e lk ^[t + l,i]{...) 
+ J2 e* w [i,* 3 ](.. J * + l I ..) 

+ X)e* fc (' +1 )[i + l > * 3 ](..,i,..) 

i 3 ^i 

+ e* Ws [i 3> i](..,* + l J ..) 

i 3 #i+l 

+ ^e lfci3 [i 3 ,i + l](..,i,..) 

i 3 ^i 

+ £ e tW » [is, <,< + !,..) • (11) 

Now we use the relation from appendix lAl 

Si • S i+ i(j,i)[i + 1, A] = + 
(Si • S i+ i) 2 (i, + = + • (12) 

By summing over all i, we obtain 

JV 

^S r S !+1 |fc) = JV |fc> + 2(oos(fe) - 1) |A> 
i=l 
JV 

^(S,-S J+1 ) 2 = iV|fc) . (13) 

i=l 

This means that for a general Hamiltonian fl} 

H\F)=E Q \F) H \k) =E \k)+2J 1 (cos(k)-l)\k) , 

(14) 

which is the well-known ferromagnetic dispersion. 

IV. MODEL AT Ji = > J 2 (6 = -§) 

One immediate consequence of the dispersion relation 
(|14[) is that the excitations become dispcrsionless at the 
phase transition point D for J\ = 

JV 

i/ = -J 2 ^(S l -S l+1 ) 2 , (15) 

i=l 

where we assume J2 < 0. At this point, it is also 
possible to express the Hamiltonian in terms of a total 
spin s = projection operator of two neighboring spins 

f()(Sj, Si_|_l) = — 3 + 3 (Si • S.; + i) 2 

JV 

H = E - 3J 2 2 P (Si, S i+1 ) , (16) 

i=l 

where i?o = —NJ2. For this model we are able to com- 
pare with known results^ which serves as an illustration 
of how to use the hws basis. We know from Eq. (fl~2"l) that 
any state with total spin s = N — 1 is an eigenstate and 



degenerate with the ground state. Excitations of the form 
\k) = exp(ikii){ii, i\ + 1}(«3, H-..i]sr) with total spin 
s — N — 2 also become eigenstates at this point. This 
can be shown again by decomposing |fc) in the following 
form 

JV 

\k) = X) etWl {ii,ii + !}(•>••,•) (17) 

ii=i 

= e M {i,i + l}(.,..,.)+e k ^{i-l,i}(.,..,.) 
+e lk ^{i + l,i + 2}(.,..,.) 
+ Yl e lH3 {i 3 ,i 3 + l}(..,i,i + l,..) . 

iaT^i— l,i,J+l 

Using Eq. (|A2|) from the appendix we find 

(Si • Si+i) 2 = 3e M {i, i + l}(. ,..,.) (18) 

+ [e lfc + e- l V"Ki + l}(.,..,.) + |fc) ■ 

Finally, we sum over all i, and obtain the dispersion re- 
lation E%{k) = — J 2 [3 + 2cos(fc)] + i?o, which is gapped 
in agreement with earlier results from Ref. pTj |. 

V. MODEL AT Jj = - J 2 < (6 = ^) 

We now turn to the other ferro-antiferromagnetic 
phase transition point E at J\ = — J 2 < 0. Using the 
basis of hws in Eq. © it is not only possible to describe 
all degenerate ferro- and antifcrromagnctic ground states, 
but also to calculate correlation functions and excited 
states exactly. The model is now given by the Hamilto- 
nian 

TV 

H = J 1 ]T(S. t .S. i+1 + (S l .S. i+1 ) 2 ) , (19) 

i=l 

for a periodic chain with even number of sites N. It 
is possible to express the Hamiltonian (|19|) in terms 
of triplet projection operator of two neighboring spins 
Pi (Si, Si+i) = (1 - Si • Si+i/2 - (Si • S l+1 ) 2 /2) 

JV 

H = S -2J 1 ^P 1 (S l ,S 4+1 ) , (20) 

i=l 

where Eq = 2NJ±. Since the operator Pi can only have 
positive eigenvalues and J\ < 0, the system has so-called 
optimum ground states^ which must be in a configu- 
ration where no two neighboring spins are in a triplet 
configuration, so that Pi|0) = for all nearest neighbors. 

Another useful representation can be obtained by us- 
ing the operator V^j+i, which exchanges the quantum 
numbers of two neighboring sites. Since both singlets and 
quintets are even under exchange, we have Pi(Si, Si + i) = 
(1 — V^j + i)/2 and therefore the Hamiltonian becomes 

JV 

H = J1N + JiJ2 V i,i+i . (21) 

i=l 
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which is in fact SU (3) invariant Hence, ground states 
of the Hamiltonian ()21[) can be constructed by requiring 
Vi,i+i |0) = |0) for all i, i.e. the ground states are in- 
variant under application of all permutation operators 
Vj^-i-ii 1 ^ 21 In order to construct such a ground state, we 
use the hws notation ([3]) to describe the antifcrromag- 
netic and the ferromagnetic ground state as follows 



\F) = (1,2,...,N-1,N) 
\AF) oc y^Qi, i2}{h, u}{h, is}- 



(22) 



where is the sum over all possible permutations of 
i'l, ...,i]\f. It is clear, that the states defined in this way 
are eigenstates of all permutations V^i+i and therefore 
immediately ground states of the Hamiltonian (|19|) with 
definite total spin s — N and s — 0, respectively. There- 
fore, the energies of the ferro- and antiferromagnetic 
states cross at the phase transition point as expected, 
but in addition we find that it is possible to construct 
ground states with any even spin s by simply combining 
permutations of quintet and singlet bonds 

|0) s oc y^(n, 12, — ,i s -i,is){is+i,i s +2}— {«jv-i,iiv} , 

(23) 

which is in accordance with SU(3) invariance. There are 
no other ground states. 

It is interesting to note that in the spin-^ chain at the 
phase transition point e in Fig. [1] the ground state can 
also be written as a permutation over all possible valance 
bonds j2£ 

For each ground state all pairs of spins are equally 
entangled with each other. Consequentially, also the cor- 
relation function is the same between any two spins in 
the chain, independent of distance 



,<0| S, • Sj |0>, = 



s(s + 1) - 2N 
N(N - 1) 



(24) 



which follows from expanding (S 2 ot ) = &i ' &j + 

'^2 i S 2 ) and can also be verified by direct calculation with 
help of the scalar product in the appendix [C] The corre- 
lation function varies from (AF\ Si ■ Si 1 ' 
to (F\ S, ■ Si 



7T=T U P 

jj \F) = 1. Higher order correlation functions 
can then also be determined itcrativcly by the use of 
P 1 (S i ,S j )\0) s = (l-S. 1 .S,/2-(S l .S J ) 2 /2)|0) s = 
for any pair of spins, so that e.g. ((Si-Sj) 2 ) =2— (Sj-Sj). 

A magnon excitation can now be constructed in a sim- 
ilar spirit as in Eq. (JTUJ) on any of the ground states. 
For example a spin-1 excitation on the antiferromagnetic 
state is given by 



I*) 



AF 



H,l2 



(25) 



Eq. (fTTj) can be used. Applying V^+i and then sum- 
ming again over all i, we obtain 

H \k) s = E x {k) \k) s E^k) = 2J 1 (cos(fc) - 1) + E , 

(26) 

where \k) is now the corresponding excitation on the 
ground state |0) s with any even spin s in Eq. (|23D . 

Interestingly, the permutational ground states in 
Eq. (|23p are in fact independent of the individual cou- 
pling strengths along the chain and are even ground 
states of higher dimensional Hamiltonians as long as 
Ji = — J2. However, the excited states in Eq. ([23]) re- 
quire translational invariance. 



VI. MODEL AT Ji = -J 2 /3 > (6 = -arctan |) 

Finally we consider the famous AKLT model^ at 
point B in the phase diagram 

N ( 1 
H = Ji ^ ( Si ■ Sj+i + -(Si ■ Sj+i) 



N 



Eq + 2Ji P2{Si, Sj+i), 



(27) 



where J 1 > 0, E = -fJVJi and P 2 (S l ,S l+1 ) = 1/3 + 
Si • Sj+i/2 + (Si • Si+i) 2 /6 is the projection operator 
onto the quintet state of two neighboring spins. It is of 
course well known how to express the ground state |0) at 
this point but it is instructive to gain an alternative 
description in terms of the hws basis. 

Analogously to the other projection points, the Hamil- 
tonian (|27jl has an optimum ground state |0), which must 
obey P 2 (Si, Si+i) |0) = for all adjacent spins. Using the 
relations in appendix IA1 we know 

P 2 (Si,S i+1 ) {i,i + l} = 

P 2 (Si,S i+1 ) ({j, i}{i + 1, k} - {j, 1 + 1}{», k}) = 

P 2 (S i ,S J+1 ) ( /A /A - - /7^\ )=0. (28) 

Accordingly, we can construct the ground state by en- 
suring that for neighboring spins which are not coupled 
by a singlet bond the corresponding crossed state is sub- 
tracted. The following state is therefore a ground state 

|0) cK^](-l) #c ™ s;imss(p) {l,2}{3,4}...{Ar-l,A^} , 

(29) 

where the number of crossings is defined as the crossing of 
singlet bonds in the graphical representation in Eq. 
For example in the N — 4 chain the ground state is given 
by 



10) 



cx ({1,2}{3,4} 




{2,4}{1,3}) 

-/7*\\ 



with k ^ 0. In order to show that this state is an eigen- 
states of the Hamiltonian (f2~l"T) . a decomposition as in 



For a finite chain with open boundary conditions, the 
Hamiltonian (|27p has a four degenerate ground states, 
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one singlet and one triplet. The singlet given by Eq. ([29]) 
with s — 0. By using the relation from the appendix 
we find for triplet states 

P 2 (S l; S l+1 ) [*,*+ 1] =0 (30) 
P 2 (S tl S*+i) + 1, k} - \j, i + *}) = . 

Hence, the triplet ground state is analogously given by 
|0) = ^(-l) #crossm9s(p) [l,2]{3,4}...{^- 1,7V} , 

(31) 

where the permutations are restricted so that the indices 
in the triplet bond always remain in ascending order. 
This ground state is a triplet of total spin s = 1 in ac- 
cordance with previous results.^ 

VII. CONCLUSIONS 

We have introduced a highest weight state basis as a 
total spin representation for the Hilbert space of spin-1 
systems. In this basis it is possible to compute scalar 
products and the action of SU (2) invariant operators for 
analytical and numerical calculations. However, the new 
basis states are not orthogonal and overcomplete. 

In the ferromagnetic phase and at the phase transition 
point D (0 = — 7r/2) it is possible to construct spin wave 
excitations explicitly in the new basis. Also the AKLT 
ground states at point B (tan 9 = —1/3) can be explic- 
itly expressed using the highest weight states. 

At the phase transition point E (0 = 37r/4) we find 
all degenerate ground states which are multiplets of even 
total spin s = 0,2, ...,N. The states can be expressed 
as permutations over all hws bonds. It is possible to 
determine the spin correlation functions (Sj-Sj) which do 
not decay along the chain, even for the antifcrromagnctic 
state. Corresponding excited states were also determined 
in the hws basis. There is some hope that these results 
at point E can be used in future works as an ansatz for 
wave- functions in the antiferromagnetic region in order 
to investigate the phase between points E and F, which 
is still not fully understood. 

Acknowledgments 

We are grateful for useful discussions with Alexander 
Struck. This work was supported in part by the Graduate 
Class of Excellence MATCOR funded by the State of 
Rheinland-Pfalz, Germany. 



APPENDIX A: ACTION OF 517(2) INVARIANT 
OPERATORS 

The action of ha+i = Sj • Sj + i on bonds where both 
indices are in the same singlet, triplet, or quintet coupling 



is immediately given by the respective eigenvalues 

hi, i+ i{i,i + l} = -2{M+1} 
h i)i+ i[i,i + l] = + 
hi, l+1 (i,i + l) = (i,i + l) . (Al) 

Using the definition in Eq. (|2|) it is straight-forward to 
determine the action of /i^i+i = Si ■ Sj+i on all possible 
hws states as follows 

l,k} = {j, i+ l}{i,k} - {j, k}{i,i+ 1} 

l,k} = \j,i + l]{i,k}-\j,k){i,i + l} 

hi,i+i(j,i){i + l,k} = + l){i,k} - (j,k){i,i + 1} 

hi,i+i(j,i)[i + l,k] = {j,i + l)[i,k} . (A2) 

The application of hf i+1 = (Si ■ S^+i) 2 is then given by 
applying S, • S,_|_i twice 

hli+i{j,i}{i + l,k} = {j,i}{i + l,k} + {j,k}{i,i + l} 
hli+i\j,i]{i + l,k} = \j,i]{i + l,k} + \j,k]{i,i + l} 

hli+i(j,i){i + l,k] = (j,i){i + l,k} + (j,k){i,i + l] 
hi, i+ i(j,i)[i + l,k} = (j,i)[i + l,k\ . (A3) 

In this way the action of any SU(2) invariant operator 
involving two spins can be calculated as a linear combi- 
nation of Si ■ S,-|_i and (Si • S^+i) 2 using above relations. 

APPENDIX B: HWS BASIS FOR N = 4 SPINS 

We consider a spin-1 system consisting of 4 sites and 
find all hws states of the form ([3]), using at most one 
triplet bond. In the spin s = sector we find 

{1,2}{3,4},{1,3}{2,4},{1,4}{2,3}; (Bl) 

in the spin s = 1 sector we find 

[1,2]{3,4},[1,3]{2,4},[1,4]{2,3}, 
[3,4]{1,2},[2,4]{1,3},[2,3]{1,4} ; (B2) 

in the spin s = 2 sector we find 

(1,2){3,4},(1,3){2,4},(1,4){2,3}, 
(3,4){1,2},(2,4){1,3},(2,3){1,4} ; (B3) 

in the spin s = 3 sector wc find 

(1,2)[3,4],(1,3)[2,4],(1,4)[2,3], 
(3,4)[1,2],(2,4)[1,3],(2,3)[1,4] ; (B4) 

and the ferromagnetic state is 

(1,2,3,4). (B5) 

Using the addition rules for angular momenta 1 ® 1 ® 
1(8)1 = 3- 6-1 6-2 © 3-3 1-4 we 
see that the states of spin s = 0, s = 1, s = 2 and 
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s = 4 are linearly independent, as can also be checked 
directly using Eq. (|2} . For the s = 3 sector we obtain the 
following linear relations 



(B6) 



APPENDIX C: SCALAR PRODUCT 

First we note that only states with equal total spin, 
i.e. states with equal number of singlet S, triplet T and 
quintet Q bonds, have a non- vanishing scalar product. 
Using Eq. ([2]) we write the hws <j3j) in the expanded form 



X t 



^(2Q + l)^I(2Q + 2) 



(ti,t 2 )6T,0i,s 2 )e£ 

,(22Q + 3) (ilQ + i) 



_i_(i2Q) 



(CI) 



,(*JV-l) 



where r = {(+, 0), -(0, +)} and E = 
{(+,-),(-,+), -(0,0)}. The algorithm how to 
compute the scalar product is based on determining the 
number of different combinations of S z product states 
in the decomposition (|Clj) that are the same in the two 
states in the scalar product. This number multiplied by 
the corresponding normalization factors gives the scalar 
product. 

For illustration let us consider two states in a chain of 
length N = 10 as an example 

|1) = (1,2,3,4)[5,6]{7,8}{9,10} 
|2) = (1,2,5,7)[3,6]{4,8}{9,10}. 

In order to calculate the scalar product (1|2) wc simply 
count the number of states that agree on both sides after 
the decomposition (|C1[) and then multiply by the nor- 
malization. The quintet spins (1, 2, 3, 4) on the left hand 
side and (1,2,5,7) on the right hand side must all be in 
the state |+). Therefore, according to Eq. ((2|) only states 
can agree on both sides of the scalar product where the 
S z eigenvalues of the spins at site 6 is |0) and at site 
8 it must be |— ). We have exactly one possibility of 
different combinations of product states to agree, which 
must be multiplied by the respective normalizations in 
Eq. ([2]). The spins with the label 9 and 10 are already 
in the same state and do not change the scalar product 



(or alternatively speaking this bond contributes a factor 
of 3 both in the numerator and denominator). There- 
fore, we obtain after multiplying by the normalization 
factors (1|2) = 3/(2 • 3 2 ) = 1/(2 • 3). Obviously, scalar 
products factorize according to " connected clusters" and 
the analysis can be done for each connected cluster sep- 
arately. From this it also follows that in order for the 
scalar product to be non-zero, at least one spin of each 
singlet bond must again be in a singlet bond on the other 
side of the scalar product. 

The application of this algorithm on connected clusters 
with s = gives a scalar product which is determined 
by the minimum number of index exchanges 7 between 
bonds in order to bring the two states into the same hws 
state 



(V>s=o(ii, -,iNMs=o(i'i, -,i' N )) = 
1 

(ip s =o{ii, *'iv)| V-V |V> s =o(«i, in)) 



(C2) 
1 

37 ' 



•} = {■■> j}{h ■} exchanges the indices 



where Vi 

between two singlet bonds. This rule can most easily 
be understood by first considering two equal states with 

(ip\ip) = 1- Then by successively changing indices with 

7 

Vij in order to create \ip') = V...V it is clear that the 
product state with 1 00) remains the same, while product 
states involving | H — ) no longer agree. Each application 
of V therefore reduces the scalar product by a factor of 
three. In a s = connected cluster of N sites, there 
are N/2 singlet bonds. The exchange operator V can at 
most be applied N/2 — 1 times between bonds in order 
to bring the states into an equivalent form. Therefore, 
the scalar product must be at least \/"i N / 2 ^ 1 or larger 
between such states. 

As an example let us consider the states 

|3) = {1,2}{3,4}{5,6}{7,8}{9,10} 
|4) = {1,10}{2,5}{4,6}{7,8}{3,9}. 

The scalar product (3|4) can be obtained from the num- 
ber of exchange operators which wc need to transform 
the state |4) into the state |3). 

(3|4) = (3|V/ 3 , 6 ^5,9^2,io|3) = ^ • (C3) 

There are other possible choices of exchange operators 
Vij, but the minimum number of exchanges is always 
7 = 3 in this case. 
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